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Abstract. Within the next few years, several instruments aiming at imaging extrasolar plan- 
ets will see first light. In parallel, low mass planets are being searched around red dwarfs 
which offer more favorable conditions, both for radial velocity detection and transit studies, 
than solar-type stars. We review recent advancements in modeling the stellar to substellar 
transition. The revised solar oxygen abundances and cloud models allow to reproduce the 
photometric and spectroscopic properties of this transition to a degree never achieved be- 
fore, but problems remain in the important M-L transition characteristic of the T c ff range of 
characterisable exoplanets. 
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1. Introduction 

Since spectroscopic observations of very low 
mass stars (late 80s), brown dwarfs (mid 90s), 
and extrasolar planets (mid 2000s) are avail- 
able, one of the most important challenges in 
modeling their atmospheres and spectroscopic 
properties lies in high temperature molecular 
opacities and cloud formation. K dwarfs show 
the onset of formation metal hydrides (start- 
ing around T eff ~ 4500 K), TiO and CO (be- 
low T e ff ~ 4000 K), while water vapor forms 
in early M dwarfs (T eff ~ 3900 - 2000 K), 
and methane, ammonia and carbon dioxide 
are detected in late-type brown dwarfs (T e g 
~ 300- 1600 K) and in extrasolar giant planets. 
Cloud formation is also an important factor in 
the detectability of biosign atures, and for the 
habitability of exoplanets dPaillet et alj 120051 
Kasting 2001). 



Send offprint requests to: F. Allard 



Extrasolar planets for which we can cur- 
rently characterize their atmospheres are either 
those observed by transit (r eff ~ 1 000- 2000 K 
depending on their radius relative to that of 
the central star) or by imaging (young plan- 
ets of r e ff ~ 500 - 2000 K depending on their 
mass and age). Several infrared integral field 
spectrographs combined with coronagraph and 
adaptive optic instruments are coming online 
before 2013 (SPHERE at the VET, the Gemini 
Planet Imager at Gemini south, Project 1640 at 
Mount Palomar, etc.). The E-ELT 41 m tele- 
scope in Spain due around 2020 will also be 
ideally suited for planet imaging. 

M dwarfs are the most numerous stars, 
constituting 80% of the stellar budget of the 
Galaxy, and around 600 brown dwarfs and 
planets are currently known despite their faint- 
ness in the solar neighborhood vicinity. Single 
very low mass (VLM) stars and brown dwarfs 
are therefore more directly observable and 
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characterizable then exoplanets. They repre- 
sent, beyond their own importance, a won- 
derful testbed for the understanding of exo- 
planetary atmospheric properties together with 
solar system studies. Planets can even share the 
atmospheric composition of brown dwarfs of 
same T e fi (see section |5]below). 

The models developed for VLMs and 
brown dwarfs are therefore a unique tool, if 
they can explain the stellar-substellar transi- 
tion, for the characterization of imaged exo- 
planets. In this paper, we review the ability of 
recently published models in reproducing con- 
straints along the M-L-T spectral transition. 



2. Model Construction 



The modeling of the atmospheres of VLMs has 
evolved (as here illustrated with the develop- 
ment of the PHOENIX atmosphere code) with 
the extension of computing capacities from an 
analytical treatment of the transfer e quation 
using moments of the radiation field 



Il960h . These are written in the compact vector 
notation as: 



quation 
(kllardl 



1990), to a line-by-l ine opacity sampl ing in 
spherical symmetry (lAllard et al.l fl994, 1997 
and Hauschildt et al. 199 9) and more recentl y 
to 3D radiation transfer (Seel mann et al.8 2010). 
In parallel to detailed radiative transfer in an 
assumed static environment, hydrodynamical 
simulations have been developed to reach a re- 
alistic representation of the granulation and its 
induced line shifts for t he sun and sun-like s tars 
(see e.g. the review bv lFrevtag et al.ll2012l) by 
using a non-grey (multi-group binning of opac- 
ities) radiative transfer with a pure blackbody 
source function (scattering is neglected). 



To illustrate the various assumptions made 
by constructing model atmospheres, let us 
begin with the description of the equations 
of ideal magnetohydrodynamics (MHD) — 
adapted here for the stellar case by specifying 
the role of gravity, radiative transfer, and en- 
ergy transport — which are themselves a spe- 
cial case (no resistivity) o f the more general 
equations (see for example lLandau & Lifshita 



= 



^ + V-(pvv + (P + ^B-B)I - BB)= pg, 

at 2 



+ V-(vB -Bv) 
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~dt 

^fl +V -((p et+ P+iB-B)v 
at 2 

-(v-B)B + F rad ) =0 . 

The vectors are noted with boldface char- 
acters, while scalars are not. For example, P 
is the gas pressure, p the mass density, g the 
gravity, and v is the gas velocity at each point 
in space. B is the magnetic field vector, where 
the units were chosen such that the magnetic 
permeability /j is equal to one. I is the identity 
matrix and a • b = 2* a kbk the scalar product 
of the two vectors a and b. The dyadic tensor 
product of two vectors a and b is the tensor 
ab = C with elements c mn = a m b„ and the nth 
component of the divergence of the tensor C is 
(V • C)„ = 2m dc mn /dx m . In this case, the total 
energy is given by 

pe t = pe-i + p^v ■ v + ^B • B + p<E> , (2) 

where e\ is again the internal energy per unit 
mass, and <t> the gravitational potential. The ad- 
ditional constraint for the absence of magnetic 
monopoles, 



V-B = 



(3) 



must also be fulfilled. 

The first, third, and last equations in eq. Q] 
correspond to the mass, magnetic field, and en- 
ergy conservation, while the second equation 
is the budget of forces acting on the gas. In the 
case of stellar astrophysics, gravitational accel- 
eration is an important source term, while the 
radiative flux participates in the energy bud- 
get. Further assumptions are made in the nu- 
merical solution of these equations to address 
different astrophysical problems in very differ- 
ent regimes. The chromospheres correspond to 
a regime of high Mach numbers and strong 
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magnetic fields where ionized gas has to fol- 
low the magnetic field lines, and where the 
radiative transfer must be solved for the case 
of a non-ideal gas. The photospheric convec- 
tion simulations correspond to a regime where 
the thermal and convective turnover timescales 
are comparable i.e. Mach numbers are around 
1, and the non-local radiative transfer must be 
solved, often for an ideal gas. And the interior 
convection and/or dynamo simulations corre- 
spond to a regime where the thermal timescale 
is much larger then the turnover timescale, 
which in turn is much larger then the acoustic 
timescale. The radiative flux can be approxi- 
mated by the diffusion approximation, and the 
magnetic field lines are dragged by the ionized 
gas. 

Radiation hydrodynamical (RHD) simula- 
tions ignore by definition the magnetic field 
terms in equation Q] This is a good approxi- 
mation when modeling the neutral photosphere 
(where most of the emitted flux emerges) of 
low mass, very low mass stars, and brown 
dwarfs — with the exception of the ultravio- 
let and visual spectral range of flaring stars and 
for the resulting emission lines. RHD simula- 
tions, especially in 2D and 3D, are computa- 
tionally expensive, and, when treating - if at all 
- radiative transfer, can currently be performed 
only for a restricted number of wavelengths, 
or wavelength bins (typically 4 to 12). In the 
case of solar-like photospheres, RHD simula- 
tions using 140 2 x 150 grid points over 5 hours 
of stellar time C05BOLD required on paral- 
lel computers (2 no des) several CPU months 
dLudwig et al.l 12009b . In the case of red and 
brown dwarf simulations, local 2D cases with a 
resolution of 400 x 30 covering 2 days of stel- 
lar time C05BOLD dFrevtag et al.l l2010h re- 
quired 1 month of CPU time. 

The classical approach for interior and 
atmosphere models consists in simplifying 
the problem for a gain of computing effi- 
ciency, neglecting the magnetic field, convec- 
tive and/or rotational motions and other multi- 
dimensional aspects of the problem, and as- 
suming that the averaged properties of stars 
can be approximated by modeling their proper- 
ties radially (uni-dimensionally) and statically. 
We also assume that the atmosphere does nei- 



ther create nor destroy the radiation emitted 
through it. Neglecting motions in modeling the 
photospheres of VLM stars, brown dwarfs, and 
planets is acceptable since the convective ve- 
locity fluctuation effects on line broadening is 
hidden by the strong van der Waals broadening 
prevailing in these atmospheres. But this is not 
the case of the impact of the velocity fields on 
the cloud formation and wind processes (see 
section [4] below). In this case, equation Q] re- 
duces to the so-called hydrostatic equation and 
constant flux approximation for the radial or z 
direction used in classical models: 

dP 

~T = -P8 > 

or 

r (4) 

dF ad = d(jF A dA) _ p 
dr dr 

This allows computing the interior evo- 
lutive properties of stars throughout the 
Hertzsprung-Russell diagram, and to solve the 
radiative transfer in the atmosphere for a much 
larger number of wavelengths (line-by-line or 
opacity sampling) or wavelength bins (Opacity 
Distribution Function or ODF, K-Coefficient) 
compared to R(M)HD simulations. Classical 
model atmospheres impose therefore the inde- 
pendent parameter F ra d (= cr Teg, where <j is 
the Stefan-Boltzman constant) and compute F^ 
so that, after model convergence, the target F ra d 
is reached. Other independent parameters are 
the surface gravity g and the abundances of the 
elements e,-. This makes it possible to create ex- 
tensive databases of synthetic spectra and pho- 
tometry that provide the basis for the interpre- 
tation of stellar observations. 

All the model atmospheres compared in 
this review are classical models in this sense, 
and differ mainly in the completeness and ac- 
curacy of their opacity database ()including 
their cloud model assumptions), and the as- 
sumed solar abundances used for the particu- 
lar grid shown. They must resolve the radiative 
transfer for the entire spectral energy distribu- 
tion (as can be seen from eq. [4]l with a good 
enough spectral resolution to account for all 
cooling and heating processes. 

Classical model atmospheres differ also 
from one another in their construction philos- 
ophy, which is linked to their period of initial 
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development. The code bv dTsui ill 19651 2 002. 
2005), the ATLAS codebv lKuruczld 19731) and 
ICastelli & Kuruczl d2004l). and the MARCS 
code bv (lGustafssonetalJll975L 2008) have 
seen the punched computer cards and the need 
to spend all efforts in saving characters and 
computer time. These models pre-tabulate their 
opacities *c,(/l) = <x,(/l), where Kt{X) 
[cm -1 ] and cr,(/l) [cm 2 ] are the opacity coef- 
ficient and cross-section at the wavelength A, 
and n, [cm -3 ] is the number density of species i 
i.e. atoms, molecules or grains) to interpolate 
them later during the mod el atmosphere exe- 
cution. The PHOENIX code dAllard et al.ll 1994 
2012) on the other hand, also to distinguish it- 
self from its forefathers, took the approach of 
computing the opacities during the model exe- 
cution (or on-the-fly). This involves computing 
the opacities for billions of atomic and molec- 
ular transitions on-the-fly, though with a selec- 
tion of the most important lines. This different 
approach makes PHOENIX much slower then 
former codes, but allows to take into account 
more consistently important physical phenom- 
ena, such as those involving a modification of 
local elemental abundances along the atmo- 
spheric structure (e.g. non-LTE, photoionisa- 
tion, diffusion and cloud formation). 

Therefore, and especially in the PHOENIX 
case, the computational requirements of clas- 
sical model atmospheres, even nowadays, pre- 
clude in practice modeling globally a star from 
its interior to its photospheric layers. Besides, 
an eventual classical static ID model appears 
less interesting then global RHD simulations. 
This is becoming possible even with rotation 
but of course at the cost of some severe ap- 
proximations at this point: the innermost core 
is replaced by an adapt ed potential function 
(Steff en&Frevtagl l2007). These RHD simula- 
tions of main sequence stars, brown dwarfs, 
and planets have also to be scaled down sig- 
nificantly in radius to preserve the ability to 
resolve convective cells and timescales of im- 
portant processes such as cloud formation. An 
alternative approach used by many authors is 
therefore to neglect small scale phenomena 
and model only larger scales, such as global 
circulation a round the planetary s urface (see 
for example iKoskinen et al.l [2007, Showman 



et al. 2009, Dobbs-Dixon et al. 2010). The 
challenge of such hydrodynamical simulations 
nowadays is to account for all the most im- 
portant opacities, in particular scattering, in 
solving the radiative transfer and hydro equa- 
tions while keeping the computing time for the 
model within reasonable limits. 

3. Molecular opacities 

While earlier work has been developed for 
the study of red giant stars, the pioneering 
work on the modeling of VLM atm o spheres 
has be en p rovid e d bv iMouldl (1 19751) . lAllardl 
(1990) and HuJ dl99ll) using a band model 
or the Just Overlapping Line Ap proximation 
(JOLA ) opacities developed by iKivel et all 
(1 19521) and ad apted for astrophysical use by 
lGoldenldl967l) . 

More realistic model atmospheres and syn- 
thetic spectra for VLMs, brown dwarfs, and 
extrasolar planets using line-by-line or opacity 
sampling techniques have been made possible 
thanks to the development of accurate opaci- 
ties calculated often ab initio for atmospheric 
layers where temperatures can reach 3000 K. 
The process of improvements was especially 
remarkable in the case of water vapor line lists. 
Indeed, water vapor has seen an important evo- 
lution through the years from band model ap- 
proximations to straight means based on hot 
flames experiments, and then to ab initio com- 
putations. Nevertheless, the atmosphere mod- 
els have failed to reproduce the strength of 
the water bands that shape the low resolution 
(R < 300) infrared spectral energy distributions 
(hereafter SED) of M dwarfs. At the lower tem- 
peratures of brown dwarfs, methane and am- 
monia rival the effect of water. The discrepan- 
cies in the model synthetic spectra were there- 
fore believed to be due to inaccurate or incom- 
plete molecular opacities. In particular water 
vapor was suspected because the discrepancies 
were observed at infrared wavelengths in the 
relative brightnesses of the flux peaks between 
water vapor bands. 

In Fig. [T]the models are compared to the 
infrared spectrum of the M8e dwarf VB10. 
One can see that the water vapor opacity pro- 
file, which shape this part of the spectrum, 



286 



Allard: modeling the stellar-substellar transition 




Wavelength (,Lim) 

Fig. 1. Fig. 1 of the review article by Allard et al. 
(2012). The near-infrared SED of VB10 is com- 
pared to synthetic spectra (T eff = 2800 K, logg = 
5.0, [M/H]= 0.0, AA = 50 A) from diverse model 
grids published through the years. All models (ex- 
cept the NextGen/UCL case) underestimate the flux 
in the K bandpass by 0.1 to 0.2 dex. 



has strongly changed over time with the im- 
provement of computational capacities and a 
better knowledge of the interaction potential 
surface. Only the most recen t ab initio results 
dPartridee & Schwenkel 19971 and the BT2 line 
list by Barber et al. 2006) confirm the ear- 
lie stJioHhimesJaboratory experiment results 
by iLudwigl (Il97ll) . Nevertheless, a lack of 
flux persist in the K bandpass in the models 
even using the most recent BT2 opacity profil e 
(e.g. the BT-Settl models ofl Allard et al.l2012l) 
Only the UCL line list dSchrvber et al] 1 19951 
due to incompleteness, and with much of its 
deviations canceling out over the bandpasses) 
could produce seemingly correct J - K colors, 
and could allow some success of this so-called 
NextGen dHauschildt et alJI 19991) model atmo- 
spheres grid in the VLM stellar regime. 

In the substellar regime, the composition 
of brown dwarfs varies rapidly with decreas- 
ing r e ff, and the variation is responsible for 
the immense change in their SED across the 
very narrow r e ff regime of the M-L-T spec- 
tral transition. If water vapor opacities only be- 



came recently reliable, this is not the case of 
the more complex methane molecule which is 
so important in brown dwarfs, and planetary 
atmospheres. The ExoMol Project supported 
by an ERC to Jonathan Tennyson (University 
College London) will allow important ad- 
vances on these fronts in the coming years. 
A new ammonia lin e list is already availa ble 
through this project (lYurchenko et al.ll201 lb . 



4. Mixing 

Stars becomes fully convective throughout 
their interior and convection reaches furthest 
out in the optically thin regions of the pho- 
tosphere in M3 and later dwarfs with r e ff be- 
low 3200 K dAllardlll990l Chabrier & Baraffe 
2000). In most model atmospheres discussed in 
this review paper, the convective energy trans- 
fer is treated us ing the Mixing Length Theor y 
(or MLT, see iKippenhahn & Weigertl 1 1994 . 
using at best a unique fixed value of the mixing 
length of 1 .0 ( 1 .25 for the ATLAS9 models, 1 .5 
for the MARCS models, etc). However, since 
convection becomes efficient in M dwarfs, the 
precise value of the mixing length matters only 
for the deep atmospheric structure and as a sur- 
fac e boundary co n dition f or inter ior models. 
Ludwigetall d2002l) and iLudwig et all 



(2006) have been able to compare the PHOENIX 
thermal structure obtained using the MLT with 
that of RHD simulations. They showed that 
the MLT could reproduce adequately (except 
for the overshoot region) the horizontally av- 
eraged thermal structure of the hydro simula- 
tions when using an adequate value of the mix- 
ing length parameter. This value has been esti- 
mated for M dwarfs to vary with surface grav- 
ity from a=llH p =\.% to 2.2 (2.5 to 3.0 for the 
photosphere). 

The BT-Settl models use the mass and 

surface grav ity dependent prescription of 

ILudwig et alJ d 19991) for hotter stars, together 
with an avera ge (2.0) of the value s derived for 
M dwarfs by (ILudwig et al]|2002l 2006). They 
use as well the micro-turbulence velocities 
from the radiation hy drodynamical simulations 
dFrevtag et al] |2010), and the velocity field 
from RH D simulations from ILudwig et al] 
(2006) and lFrevtagetal.l(l2010h to calibrate the 
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scale height of overshoot, which becomes im- 
portant in forming thick clouds in L dwarfs but 
is negligible for the SED of VLMs and brown 

dw arfs otherwise . 

iFrevtag et~ai1(l2010l) have indeed addressed 
the issue of mixing and diffusion in VLM 
atmospheres by 2D RHD simulations, using 
the PHOENIX gas opacities in a multi-group 
opacity scheme, and forsterite with geometric 
cross-sections. These simulations assume effi- 
cient nucleation, using initial monomer densi- 
ties estimated from the total available density 
of silicon (least abundant element in the so- 
lar composition involved in forsterite). They 
found that gravity waves form at the internal 
convective-radiative boundary, and play a de- 
cisive role in cloud formation, while around 
T e ff « 2200 K the cloud layers become thick 
enough to initiate cloud convection, which 
dominates in the mixing. 

5. Atmospheric composition 

The composition of the atmospheres of stars, 
brown dwarfs, and planets is a function of T e ff 
(radiation either due to internal heat from nu- 
clear fusion and contraction or from irradiation 
by a parent star), of surface gravity to a lesser 
extent, and of the elemental abundances of the 
initial gas from which the star or stellar system 
is formed. Stellar model atmospheres assume 
scaled solar abundances for all elements rel- 
ative to hydrogen. Additionally, some enrich- 
ment of a-process elements (C, O, Ne, Mg, Si, 
S, Ar, Ca, and Ti) resulting from a "pollution" 
of the star-forming gas by the explosion of a 
supernova is appropriate in the case of metal- 
poor subdwarfs of the Galactic thick disk, halo, 
and globular clusters, and the stars in the high 
stellar den sity environment to wards the galac- 
tic center dGaidos et al.ll2009l) . 

5.1. The revision of solar abundances 

Important revisions have been made to the so- 
lar abundances based on radiation hydrody- 
namical simulations of the solar photosphere, 
and to improvements in the detailed line pro- 
file analysis. Indeed, two separate groups us- 
ing independent RHD and spectral synthesis 



codes (lAsnlund et al.ll2009l Caffau et al. 201 1) 
obtain an oxygen reduction of 0. 1 1 - 0. 19 dex 
(up to 34%) c ompared to the previ ously used 
abundances o flGrevesse et al.l(ll993l) . Since the 
overall SED of late K dwarfs, M dwarfs, brown 
dwarfs, and exoplanets is governed by oxy- 
gen compounds (TiO, VO in the optical and 
water vapor and CO in the infrared), the el- 
emental oxygen abundance is of major im- 
portance. Fig. [2] shows an example of these 
effects, where several models are compared 
to the optical to infrared SED of the M5.5, 
M9.5, and L0 dwarfs of the LHS 1070 sys- 
tem. The BT-Setfl model by Allar d et al (2012) 
is based on the lAsplund et alj d2009f) solar 
abundance v a lues, w hile DRI FT models by 
iHelling et ail (l2008bl) use the IGrevesse et all 
d 1993b s olar abundances and the MARCS 
model by Gustafsson et al. (2008) uses the val- 
ues o flGrevesse et al.1 (120071) . Inspecting Fig. [2] 
one can see that the MARCS model show a 
systematic near-infrared flux excess, compared 
both to observations and the other models, 
which is probably caused by t he much lower 
oxyge n abundance values of IGrevesse et al.1 
(2007). The oxygen abundances sensitivity of 
TiO bands is expressed as a reduced line 
blanketing effect at longer wavelengths, par- 
ticipating in the w ater vapor profile changes 
(lAllardetal.1120001) . 

The influence of the solar oxygen abun- 
dance can als o be clearly seen i n Fi g. [3] which 
compares the lCasagrande et alj (120081) T P ff and 
metall icity estimates with the iBaraffe et al.l 
(1998) NextGen isochrones (assuming an age 
of 5 Gyrs) using model atmospheres from var- 
ious authors. The oxygen abundance effects 
are particularly highlighted by comparing the 
BT-Se tfl model based on the Asplu nd et alj 
(2009) values with models based on earlier so- 
lar abundance values. This is the case of the 
AMES-Cond/Dusty and BT-NextGen models 
by Allard et al. (2001, 2012) which are based 
on the Grevesse et al. (1993) solar abundances. 
On can see that the higher oxygen abundance 
causes models to appear too blue by as much 
as 0.75 mag compared to models based on 
the I Asplund et all ( 120091) values. The MARCS 



models ( Gustafsson et al. 2008) based on the 
IGrevesse et al.1 ( 20071) values show on the con- 
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Wavelength (A) Wavelength (|im) 

Fig. 2. Fig. 8 (left) and 9 (right) of Rajpurohit et al. (2012). In each of the three panels, the spectra of the 
resolved components of the L HS 1070 triple system dLeinert et ali l994) — observed with the Faint Object 
Spectrograph (FOS) on the left dLeinert et al.l2 000) and with NICMOS on HST on the right — are compared 
to models of various authors. Black: observed spec tra. Green: best^ 2 fit BT-Settl AGSS model Allard et a l. 
(2012). Blue: DRIFT model l lHelling et al]|2008bl). And brown : MARCS model dGustafsson etaill2008l) . 
The obtained parameters are given in iRajpurohit et al.l d2012l) . but do not change by more then one to 
two si gma from earlier estimations (Leinert et al. 2000J) based on the AMES-Dusty models bv I Allard et alj 
(2001), given the grid parameter spacing (100 K). The B and C components are cool enough (r eff = 2500 
and 2400 K respectively) for their SED to be affected by dust formation. 



trary a systematically increasing excess in J - 
K s with decreasing T e ff. The models are most 
sensitive on the solar oxygen abundances for M 
dwarfs around 3300 K, i.e. at the onset of water 
vapor formation. 



Th e NextGen model by lHauschild t et alj 
d 19991) dates too far back and suffers from too 
much opacity differences (incompleteness es- 
sentially) to participate in this illustration. In 
fact, this plot helps to conclude that using the 
NextGen models caused a systematic overesti- 
mation of for VLM stars. It is interesting 
to note that all models appear too red in the K 
dwarf range above 4000 K. This may be due 
to an under representation of the K dwarfs in 
this diagram. The unified cloud model (here- 
after UCM) bv lTsuiil d2002l) show a completely 
different behavior in this diagram, sharing the 
colors of NextGen or even MARCS models at 
4000 K, but diverging towards the BT-Settl col- 
ors at 3500 K to finally cross-over to bluer col- 
ors as dust begin to form and affect the SED 
below 2600 K. 



The various model atmospheres have not 
been used as surface boundary condition to 
interior and evolution calculations, and sim- 
ply provide the synthetic color tables interpo- 
lated on the publish ed theoretical isochrones 
dBaraffeetal.lll998l) . Even if the atmospheres 
partly c ontrol the cooling and evo lution of M 
dwarfs (Chabr ier & Baraffdl 19971) . differences 
introduced in the surface boundary conditions 
by changes in the model atmosphere composi- 
tion have negligible effect. 



5.2. Metallicity 

The ground work in understanding the metal- 
licity effects on the SED and colors o f VLM 
star s has been establish e d bv lAllardl (Il990l) 
and lAllard & Hauschildtl d 19951) . and summa- 
rized bv lAllard et alj (Il997l) . The main effects 
of reducing the metallicity are the gradual dis- 
appearance of the double-metal molecules of 
importance for the overall opacities and the 
SED of VLM stars (TiO, VO, CO), and the 
increased pressure effects (atomic line widths, 
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Fig. 3. Estimated T eB and metallicity (lighter to darker tones) for M dwarfs by Casagra nde et al.l (2008) 
on the left, and brown dwarfs by Golimowski et al. (2004) and V rba et al. (2004) on the right are com- 
pared to the NextGe n isochrones for 5 Gyr s Igaraffe et al. (1998) using model atm ospheres by various 
aut hors: MARCS bvlGustafsson et al] d2008h. AT LAS9 bv lCastelli & Ku ruczl J2004I) DRI FT-PHOENIX 
by iHelling et al] <2008bl) . UCM by iTsuiil i200$>. Clear/Cloudy b y iBurrows et alN2006l) . NextGen by 
iHauschildt et al] dl999h . AMES-Cond/Dustv by lAUard et al.l d200lh . and the BT models by Allard et al. 
(2012). The region below 2900 K is dominated by dust formation. The dust free models occupy the blue 
part of the diagram and only at best explain T dwarf colors, while the Dusty and DRIFT models explain at 
best L dwarfs, becoming only redder with decreasing r eff . The BT-Settl, Cloudy and UCM r crit = 1700 K 
models describe a complete transition to the red in the L dwarf regime before turning to the blue into the T 
dwarf regime. The Cloudy model however does not explain the reddest L dwarfs. 



the strength of the hydride molecular bands in- 
cluding the well-known H2 collision-induced 
absorption (CIA) bands in the K bandpass) re- 
sulting from the increased transparency of the 
atmosphere. This is illustrated by Fig. 5 of 
Alla rd et al] (1 19971) . which shows how these 
changes operate in pushing the peak of the 
SED to the blue as metallicity decreases from 
solar to [M/H]=-2.0, and -4.0, and the optical 
is becoming brighter while the CIA opacities 
depress the K-bwA flux. 

The main difficulty in parametrizing M 
subdwarfs using pressure indicators are due to 
the fact that an atmospheric pressure increase, 
while obtained by decreasing the metallicity, 
can also be provided by increasing the surface 
gravity, and/or reducing the effective tempera- 
ture. Moreover, subdwarfs have systematically 
higher gravities and smaller radii than solar 
type M dwarfs (Bar affe et al.lll997l) . compen- 
sating for the increased brightness of the more 
transparent atmosphere at optical wavelengths. 
Thus, disentangling these effects can be diffi- 



cult to impossible using low spectral resolution 
or broad band colors, and detailed high resolu- 
tion studies are necessary. Work is in progress 
using the zeta index (or TiO/CaH flux ratio) by 
Lepine & Scholz (2008) and SDSS g-r and r- 
i colors (Barbara Royas Ayala and Sebastien 
Lepine, private comm.). 

In the substellar regime, the metallicity ef- 
fects are similar to those enumerated above, 
but adapted to their cooler composition, with 
the difference that the SED of T dwarfs is 
not expected to become bluer with decreas- 
ing metallicity as for M dwarfs. This is due to 
the already extreme transparency of their atmo- 
sphere (due to condensation), and the result- 
ing extreme strength and width of optical alkali 
spectral doublets (in particular Na I D and K I 
at 0.77/im). 

Even though some metal-poor L dwarfs 
have been identi fied by Burgasser et a l. (2003, 
2004, 2006) and lSivaranietall d2009l) . metal- 
poor brown dwarfs have all chances to be T 
dwarfs by the age of the thick galactic disk 
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(Population II) according to evolutionary mod- 
els. For this reason, and because T dwarfs are 
more readily modeled than the dustier L dwarfs 
(see section [6] below) , severa l studies such 
as tho se bv lLiuetal.1 ( 2007). Legge ttet al.1 
d2010l) . and iPinfield et alj d2012l) have been 
done, despite the extreme difficulty tied to the 
poor spectral resolution of the observations and 
the fact that the SED of T dwarfs changes very 
slightly with the atmospheric parameters com- 
pared to those of L and M dwarfs. The sub- 
metallicity of L and T brown dwarfs is not 
extreme ([Fe/H] > -0.5). More metal-poor 
brown dwarfs would more likely be Y dwarfs. 



5.3. Disequilibrium chemistry 

Mixing has been held responsible, beyond its 
role in updrafting condensible gas into the 
cloud forming layers above, for deviations 
from gas phase chemical equilibrium in the 
atmospheres of ultracool T dwarfs, exoplan- 
ets, and Jupiter. These effects on their SED 
have been observed as an excess of carbon 
monoxide absoiption (iNoll et al.ll 19971 Griffith 
& Yelle 1999). Similarly, amm onia has been 
show n to be under-abundant (Saum on et al.l 
2006). More recently, carbon dioxide, which 
under chemical equilibrium conditions is ex- 
pected to form in only very small quantities 
in hydrogen-rich atmospheres, ha s been di- 
rectly detected spectroscopically (Ts uii et all 
1201 1|) and inferred from mid- infrared photom- 
etrv dBurningham et al.ll20T"fl) . 

This is understood as the result of the fact 
that the formation of methane and ammonia, 
which are expected to dominate the carbon 
and nitrogen chemistries, respectively, in the 
low te mperature limit under eq uilibrium condi- 
tions (jLodders & Feglevl2006l) . has to compete 
against upmixing from the deeper and warmer 
atmospheric layers. If the local mixing rate is 
high compared with the relevant reaction rates, 
their high-temperature counterparts (CO, CO2, 
N2) can instead be observed in the upper atmo- 
sphere in excess of their local chemical equi- 
librium concentrations. To estimate the forma- 
tion timescales one needs to identify the most 
efficient reaction path of formation and then 
isolate the rate-limiting, i. e. slowest, step in 



this path. For the conversion of nitrogen to am- 
monia this net reaction and limiting step is gen- 
erally taken to be 



3N 2 + 2H 2 = NH 3 : 
N 2 + H 2 ^ 2NH 



(5) 



according to Lewis &Prinnl (Il980h . with the 
resulting timescale given by 



f^ 1 = 8.45 x 10-V 8151/7 -™ 3 



cmV [H 2 ] 



(6) 



where T is the gas temperature in K and [H 2 ] 
the number density of molecular hydrogen. 
Due to the strong temperature sensitivity this 
makes nitrogen destruction quickly inefficient 
around 2000 K, meaning that the nitrogen- 
ammonia ratio in typical brown dwarf atmo- 
spheres is fixed already in the deep convection 
zone. 

The case for the carbon monoxide to 
methane conversion is more complex, where 
a variety of possible reaction paths and corre- 
sponding timescales have been discussed in the 
literature. In an extensive analysis of the reac- 
tion network fVisscher et al.l (1201 Ol) propose 



CO + 3H 2 = CH 4 + H 2 : 

H 2 + CH3O ^ CH3OH + H 

with the timescale given by 



l co 



*R863 - 



[H 2 ][CH 3 Q] 
[CO] 



(7) 



(8) 



and estimate a reaction constant 
£r863 = 1.77xlO- 22 r- 3 09 e- 3055/:r cm 3 s _1 .The 
resulting timesca les are at some varianc e with 
earlier estimates dPrinn & Barshavll 19771 Yung 
et al. 1988, Griffith & Yelle 1999, and Schaefer 
& Fegley 2010), but they are generally becom- 
ing longer than atmospheric mixing timescales 
only in the cooler outer layers above the con- 
vection zone. This means that measurements 
of the CO abundance are a sensitive probe of 
overshoot and other mixing processes above 
the Schwarzschild boundary. 

This mixing may be described by a sin- 
gle diffusion coefficient introduced as an ad- 
ditional model parameter that can be in- 
ferred from the height above which CO 
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and CH4 should be kept fixed at their rela- 
tive abundances; however this method strictly 
only probes mixing at this "quench level" 
(ISaumon & Marlevll2008l). T he RHD simula- 
tions of iFrevtag et al.1 (1201 ol) have allowed to 
model the underlying mixing processes as a 
function of height, although the translation 
from the averaged hydrodynamic velocity field 
to the molecular diffusion coefficient is still 
subject to some uncertainty: the scaling be- 
havior of stochastic waves (in stable regions 
of the atmosphere, far away from the con- 
vection zones) is different from that of de- 
veloped turbulence (as found - approximately 
- in the deep stellar convective layers or the 
thin cloud convection zone) or from motions 
in ove rshoot layers. Fig. 13 of Frevtag et al. 
(120101) shows how the former may be approxi- 
mated by scaling with some power of the Mach 
number of the flow. Spectroscopic comparison 
with the strength of CO and CH4 fea tures, e. g. 
in Figs. 25, 26 of lKing et al] ( 1201 Ol) demon- 
strates that these latter chemical models with 
relatively fast rates are most consistent with 
observations. 



6. Cloud formation 

One of the most important challenges in 
modeling t hese atmospheres i s the formation 
of clouds. iTsuii et all (1 19961) had identified 
dust formation by recognizing the condensa- 
tion temperatures of hot dust grains (enstatite, 
forsterite, corundum: MgSi03, Mg2Si04, and 
AI2O3 crystals) to occur in the line-forming 
layers (t * 10~ 4 - 10~ 2 ) of their models. 
The onset of this phase transition occurs in M 
dwarfs below T e ff = 3000 K, but the cloud lay- 
ers are too sparse and optically thin to affect 
the SED above T e ff = 2600 K. The cloud com- 
position, according to equilibrium chemistry, is 
going from zirconium oxide (Zr02), refractory 
ceramics (perovskite and corundum; CaTiO^, 
AI2O3), silicates (e.g. forsterite; Mg2Si04), to 
salts (CsCl, RbCl, NaCl), and finally to ices 
(H 2 0, NH 3 , NH 4 SH) as brown dwarfs cool 
down over tim e from M through L , T, and Y 
spectral types d Allard et al.l 1200 ll Fergley & 
Lodders 2006). This crystal formation causes 
the weakening and vanishing of TiO and VO 



molecular bands (via CaTi03, TKD2, and VO2 
grains) from the optical spectra of late M 
and L dwarfs, revealing CrH and FeH bands 
otherwise hidden by the molecular pseudo- 
continuum, and the resonance doublets of al- 
kali transitions which are only condensing onto 
salts in late-T dwarfs. The scattering effects 
of this fine dust is Rayleigh scattering which 
provides veiling to the optical SED, while the 
greenhouse effect due to the dust cloud causes 
their infrared colors to become extremely red 
compared to those of hotter dwarfs. The up- 
per atmosphere, above the cloud layers, is de- 
pleted from condensible material and signifi- 
cantly cooled down by the reduced or missing 
pseudo-continuum opacities. 

One common approach has been to explore 
the limiting properties of cloud formation. One 
limit is the case where sedimentation or gravi- 
tational settling is assumed to be fully effi cient. 
This i s the case of the Case B model of Tsuii 
d2002l) . the AMES-Cond mo del of Allard et 



al. (20 01), the Clear model of Burgasser et al 
(2002), and the Clear model of Burrows et al 



(2006). The other limit is the case where gravi- 
tational settling is assumed inefficient and dust, 
often only forsterite, forms in equilibrium with 
the gas p hase. This is the case of of the Case A 
model oflTsuiil d2002l), th e AMES -Dusty mod- 
els of lAllard et al.1 (12001). the BT-Dusty mod- 
els of Allard et al. (20 12), the Dusty model of 



Burgasser et al. 



2002), and the Cloudy model 



of Burrows et al. (2006). To these two limiting 
cases we can add a third case also explored by 
several, which is the case where condensation 
is not efficient and the phase transition does 
not take pl ace. This is the c a se of t he NextGen 
models of lHauschildt et alJ dl999l) . of the BT- 
NextGen models of Allard et al. (2012), and 
the Case B models of (Tsuii 2002, not shown). 

The purpose of a cloud model is to go be- 
yond these limiting cases and define the num- 
ber density and size distribution of condensates 
as a function of depth in the atmosphere, and 
as a function of the atmospheric parameters. 
The discovery of dust clouds in M dwarfs and 
brown dwarfs has therefore triggered the de- 
velopment of cloud models building up on pio- 
neering work in the co ntext of planet a ry atmo- 
spheres developed by iLewisI d 19691) . iRossowl 
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dl978l) . and iLunine et alJ dl989h . The Lewis 
model is an updraft model (considering that 
condensation occurs in a gas bubbles advected 
from deeper layers). By lack of knowledge of 
the velocity field and diffusion coefficient of 
condensates in the atmospheres of the planets 
of the solar system, Lewis simply assumed that 
the advection velocity is equal to the sedimen- 
tation velocity, thereby preserving condensible 
material in the condensation layers. This cloud 
model did not account for varying grain sizes 
(these naturally vary as a function of depth in 
the cloud layers). Rossow, on the other hand, 
developed characteristic timescales as a func- 
tion of particle size for the main microphysical 
processes of importance (condensation, coag- 
ulation, coalescence, and sedimentation). The 
curve intersections gives an estimate of the 
condensate number densities and mean grain 
sizes. However, this model made several ex- 
plicit assumptions concerning the efficiency of 
supersaturation, the coagulation, etc. 

iHeiling et alJ d2008al) have compared dif- 
ferent cloud models and their impact on model 
atmospheres of M and brown dwarfs. Most 
cloud models define the cloud base as the evap- 
oration layer provided by equilibrium chem- 
istry. In the unified cloud model of Tsuji et 
al. (2002, 2004) a parametrization of the ra- 
dial location of the cloud top by way of an ad- 
justable parameter T cr \ t was used. This choice 
permits to parametrize the cloud extension ef- 
fects on the spectra of these objects without re- 
solving the cloud model equations. In princi- 
ple, this approach does not allow to reproduce 
the stellar-substellar transition with a unique 
value of r cr it since the cloud extension depends 
on r e ff. Indeed, the transparent T dwarf atmo- 
spheres can only exist if the forsterite cloud 
layers retract below the line-forming regions in 
those atmospheres. 

lAckerman & Marlevl (|2001|) have solved 
the particle diffusion problem of condensates 
assuming a parametrized sedimentation effi- 
ciency / se d (constant through the atmosphere) 
and a mixing assumed constant and fixed to 
its maximum val ue (maximum of the inner 
convection zone ). Burgasser et al. (2002) and 
ISaumon & Marlevl (120081) found that their so- 
called Cloudy models could not produce the 



M-L-T spectral transition with a single value 
of /sed- This conclusion prompt ed them to 
propos e a patchy cloud model Marle veit al.l 
(2010). We have not been able to obtain these 
models for comparison in this paper. 

lAllard et al.l (120031) and Allard et al. (2012) 
have developed PHOENIX version 15.05 using 
the index of refraction of 55 condensible 
species, and a slightly modified version of the 
Rossow cloud model obtained by ignoring the 
coalescence and coagulation, and computing 
the supersaturation consistently. They density 
and grain size distribution with depth in the 
atmosphere is obtained by comparing the 
timescales for nucleation, condensation, gravi- 
tational settling or sedimentation, and mixing 
derived from the Mixing Length Theory 
for the convective mixing in the convection 
zones, exponential over shoot according to 



iLudwig et ai] d2002L 120061). and from g ravity 
waves according to Freytag et alJ (l2010h . The 



cloud model is solved layer by layer inside 
out (bottom's up) to account for the sequence 
of grain species formation as a function of 
cooling of the gas. Among the most important 
species forming in the BT-Settl model are 
Zr0 2 , A1 2 3 , CaTiC-3, Ca 2 Al 2 Si0 7 , MgAl 2 4 , 
Ti 2 3 , Ti 4 7 , Ca 2 MgSi 2 7 , CaMgSi 2 6 , 
CaSi0 3 , Fe, Mg 2 Si0 4 , MgSi0 3 , Ca 2 Si0 4 , 
MgTiC-3, MgTi 2 5 , Al 2 Si 2 Oi3, VO, V 2 3 , 
and Ni. At each step, the gas phase is adjusted 
for the depletion caused by grain formation 
and sedimentation. The grain sizes (a unique 
maximum value per atmospheric layer) are 
determined by the comparison of the different 
timescales and thus varies with depth to reach 
a few times the interstellar values (used in the 
dusty limiting case models) at the cloud base 
for the effective temperatures discussed in this 
paper. While the BT-Settl model assumes dirty 
spherical grains in the timescales equations to 
calculate the growth and settling of the grains, 
it only sums the opacity contributions of each 
species in each layer as for an ensemble of 
pure spherical grains. 

iHelling et alJ d2008bh and IWitte et alJ 
(2009) modified the PHOENIX code to compute 
the DRIFT-PHOENIX models, considering the 
nucleation of only seven of the most important 
solids (Ti0 2 , A1 2 3 , Fe, SiQ 2 , MgO, MgSi0 3 , 
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Mg2SiC>4) made of six different elements. The 
cloud model is based on resolving the moment 
equations for the dust density accounting 
for nucleation on seed particles and their 
subsequent growth or evaporation, solving 
from top to bottom of the atmosphere. This 
model assumes dirty grains mixed according 
to the composition of each atmospheric layer. 
It uses composite optical constants resulting 
in absorption and scattering properties of the 
grains that are therefore different than those of 
the BT-Settl models, possibly producing more 
opaque clouds. However, since the opacities 
are dominated by atomic and molecular 
opacities over most of the spectral distribution 
in this spectral type range, the impact of those 
differences are difficult to identify. The largest 
differences between the BT-Dusty, BT-Settl 
and DRIFT models are the differences in the 
local number density, the size of dust grains, 
as well as their mean composition, which are 
the direct results of the cloud model approach. 
The DRIFT model includes, similarly to the 
BT-Settl model, mixing by convection and 
overshooting by assuming an exponential 
decrease in mass exchange frequency in the 
radiative zone. But it neglects the contribution 
of the gravity waves included in the BT-Settl 
model. 

The models using the limiting cases of 
maximum dust content describe adequately 
(given the prevailing uncertainties) the infrared 
colors of L dwarfs. The cloud-free limiting 
case models, on the other hand, allow to repro- 
duce to some degree the colors of T dwarfs. 
But pure equilibrium chemistry models with- 
out parametrization of the cloud extension in 
the atmosphere cannot reproduce the observed 
behaviour of the M-L-T transition, the dusty 
models only becoming redder and dustier with 
decreasing T e ff, while dust-free models miss 
completely the reddening due to the dust 
greenhouse effects in the L dwarf regime. 
Fig. [3] shows this situation compared with the 
effective temperatures estimates obtained by 
integration of the observed SED (Golimowski 
et al. 2004, Vrba et al 2004). One can see from 
Fig. [3] that the late-type M and early-type L 
dwarfs behave as if dust is formed nearly in 
equilibrium with the gas phase with extremely 



red colors in some agreement with the AMES- 
Dusty models. The BT-Settl models (full black 
line) reproduce the main sequence down to the 
L-type brown dwarf regime, before turning to 
the blue in the late-L and T dwarf regime as 
a result of the onset of methane formation in 
the K s bandpass. The BT-Settl models succeed 
as good as the limiting case AMES-Dusty (full 
red curve), BT-Dusty (dashed red curve), and 
UCM TVrit = 1700K (full red with big dots 
curve) at explaining the reddest colors of L 
dwarfs (assuming an age of 5 Gyrs). The fact 
that a UCM model with T ciit value of 1700K 
succeeds rather well in reproducing the L-T 
transition suggest that the cloud extension is 
somewhat constant through that transition. The 
DRIFT models, on the other hand, (magenta 
with diamonds curve) reach slightly less to the 
red and do not extend low enough in temper- 
ature to explain the L-T transition. The M-L 
transition is not reproduced by any of the dif- 
ferent models, as shown by Fig. [2] where the 
CIFIST and BT-Settl models begin to show a 
/-band flux excess for the B and C compo- 
nents. This suggests that an additional element 
neglected thus far is at play, such as larger 
maybe porous grains. Indeed, all models as- 
sume thus far spherical and non-porous grains. 
The choice of solar abundances and the com- 
pleteness of the opacity databases used is also 
somewhat important. One sees in Fig.[3] (on the 
right) that models based on the Asplun det al.1 
(2009) solar abundances reach to redder col- 
ors in better agreement with constraints above 
2000 K then other models. The understanding 
of the M-L transition between r e ff = 2000 and 
2400 K is an extremely important regime for 
the study of extrasolar planets. . . 

7. Conclusions 

We have compared the behavior of the recently 
published model atmospheres from various au- 
thors across the M-L-T spectral transition from 
M dwarfs through L type and T type brown 
dwarfs and confronted them to constraints. If 
the onset of dust formation is occurring be- 
low T e s - 2900 K, the greenhouse or line blan- 
keting effects of dust cloud formation impact 
strongly (J - K s < 2.0) the near-infrared SED 
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of late-M and L-type atmospheres with 1300 < 
r eff < 2600 K. The BT-Settl models by Allard 
et al. (2012) are the only models to span the 
entire regime. In the M dwarf range, the re- 
sults appear to favor the BT-Settl based on the 
lAsplund et al.l (2009) solar abundances versus 
MARCS and ATLAS 9 models based on other 
values. In the brown dwarf (and planetary) 
regime, on the other hand, the unified cloud 
model by Tsuii ( 2002) succeeds extremely well 
in reproducing the constraints, while the BT- 
Settl models also show a plausible transition. 
However, no models succeed in reproducing 
the M-L transition between 2900 and 2000 K. 
This r e ff range is similar to that of young (di- 
rectly observable by imaging) and strongly ir- 
radiated planets (Hot Jupiters). 
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